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Abstract 

The Fokkcr-Planck equations describe time evolution of probabil- 
ity densities of stochastic dynamical systems and are thus widely used 
to quantify random phenomena such as uncertainty propagation. For 
dynamical systems driven by non-Gaussian Levy processes, however, 
it is difficult to obtain explicit forms of Fokker-Planck equations be- 
cause the adjoint operators of the associated infinitesimal generators 
usually do not have exact formulation. In the present paper, Fokker- 
Planck equations are derived in terms of infinite series for nonlinear 
stochastic differential equations with non-Gaussian Levy processes. A 
few examples are presented to illustrate the method. 

PACS 2010: 05.40.Ca; 05.10.Gg; 02.50.Fz 
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1 Introduction and statement of the problem 

The Fokker-Planck equations are one of the deterministic tools to quantify 
how randomness propagates or evolves in nonlinear dynamical systems. For 
stochastic differential equations (SDEs) with Gaussian processes such as 
Brownian motion, the Fokker-Planck equations are well established [6j[3]. 
However, for SDEs with non-Gaussian processes such as Levy processes, 
explicit forms of the Fokker-Planck equations are not easily available, except 
in some special cases [U [9]. The difficulty is to obtain the expressions for 
the adjoint operators of the infinitesimal generators associated with these 
SDEs. 

Levy processes are a class of stochastic processes having independent 
and stationary increments, as well as stochastically continuous sample paths 
[TJ [8]. Given a sample space $7, together with a probability measure P and 
the corresponding mathematical expectation E. A Levy process Lt, taking 
values in R d , is characterized by a drift parameter b £ M. d , a d x d positive- 
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definite covariance matrix A and a measure v defined on M, d and concentrated 
on M d \{0}. In fact, this measure v satisfies the following condition [T] 

|2/|| 2 A l)v{dy) < oo, (1) 

JR d \{0} 

or equivalently 

|2 



\y 
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i^(dy) < oo. 



l + lll/l 

Here || • || is the usual Euclidean norm or length in K d . This measure v is 
called a Levy jump measure for the Levy process Lf. A Levy process with 
the generating triplet (b, A, v) has the Levy-Ito decomposition 

L t = bt + B t + f yN(t,dy)+ I yN(t,dy), (3) 
J \\y\\<i J \\y\\>i 

where N(dt, dx) is the Poisson random measure, N(dt, dx) = N(dt, dx) — 
v{dx)dt is the compensated Poisson random measure, and Bt is an indepen- 
dent (i-dimensional Brownian motion (i.e., Wiener process) with covariance 
matrix A. Equation ([3|) can be formally rewritten in a differential form as 

dL t = bdt + dB t + ( yN(dt, dy) + ( yN(dt, dy). (4) 
■^IMI<i J \\y\\>i 

We shall consider stochastic dynamical systems described by the follow- 
ing SDE in the ltd form 

dX t = f(X t ,t)dt + a(X t „,t)dL t , X = x, (5) 

or in the Marcus form 

dX t = f(X t ,t)dt + a(X t _,t)odL t , X = x, (6) 

where Lt is a Levy process with the generating triplet (6, A, v). Equation 
([5]) is meaningful in the sense of 

rt rt 



X t = X + [ f(X s ,s)ds+ [ a(X s _,s)dL s , 
Jo Jo 



(7) 
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where the last term is an Ito integral. Note that Levy processes are semi- 
martingales, and Ito integrals with respect to semimartingales are well de- 
fined [7]. The equation ([6]) is interpreted as 

X t = X + I f(X s ,s)ds+ [ a(X s .,s)odL s , (8) 
Jo Jo 

where "o" indicates Marcus integral [U Q] defined by 

f <x(X s -,s)odL s = f a(X s -,s)dL s + l f a(X s _,s)a'(X s ^s)d[L s ,L s } c 
Jo Jo * Jo 

+ K(AL(s),a(X s _, S ),X s _)-X s _-a(X s _, S )A J L s ], 

0<s<t 

(9) 

with £(r,g(x),x) being the value at z = 1 of the solution of the following 
ordinary differential equation: 

^-y(z)=rg(y(z)), y(0) = x. (10) 

Assumption (Hi): 

We assume that the appropriate Lipschitz and growth conditions on the 
drift / and noise intensity a are satisfied so that both equations ([5]) and ([6]) 
have unique solutions. 

The main objective of this paper is to derive an expression of the Fokker- 
Planck equations for nonlinear dynamical systems described by SDEs (J5]) 
and ©, in ^}2] and £j3j respectively. For simplicity, we only consider one- 
dimensional case, i.e., Xt and Lt in ([5]) and ([6]) are all scalar processes. The 
conclusion can be generalized into higher dimensional cases which, however, 
will not be considered in this paper. A few examples will be presented in 

SI 

2 Fokker-Planck equations for SDEs with Ito in- 
tegrals 

In this section, we derive the Fokker-Planck equation for Ito SDE ([5|), in 
which L t being a scalar Levy process with the triplet (b, A, v). In this case, 
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A is a non- negative scalar. 

Substituting dU into (0), we get 



dX t = f(X t ,t)dt + ba(X t .,t)dt + a(X t _,t)dB t + [ a(X t -,t)yN(dt, dy) 

Av\<i 

a(X t _,t)yN{dt,dy). 



+ 
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l>i 



By ltd formula [T], it follows from (jlip that for any smooth function 4>(x), 
<P(X t+At ) - 4>(X t ) = f(X e -,s)-^-4>(X s )ds + J t+At ba(X s -,s)-^cj>(X s )ds 

rt+At ft a rt+At p,2 

/ a (X s .,s) — 4>(X s )dB s + - j a 2 (X s .,s) — ^(X s )ds 

rt+At ,■ 

/ / [<f,(X.-+ V a(X.-,a))-4>(X.-)] N(ds,dy) 
Jt J |y | > 1 

[<j>{X s . + y*(X^,s))-<t>(X a _)) N(ds,dy) 



+ 



+ 



+ 



+ 



t+At 



t J\y\<i 
rt+At r 

Jt J\y\<l 



<t> (X s _ + ya(X s .,s)) - 0(X S _) - ya(X s -, a) — 0(X S _) 

(12) 



^(dj/)ds . 



Let Q be the infinitesimal generator associated with the Markovian solu- 
tion process Xt- Then for any function cb in C^°(]R), consisting of smooth 
functions with compact support on M, we obtain 
E{0(X t+At )\X t = x)}-4>{x) 



Qd>(x) = lim 
y ' At-+0 



At 



d d A d 2 



|y|<i 



lwl>i 



d 

>0 + ya{x, t)) - <j)(x) - ya(x, t)—c/)(x) 
(x + ya(x,t)) -4>(x)] v(dy), 



u(dy) 



(13) 



or equivalently, 

d d A d 2 

Q(j){x) = f(x,t)—(f>(x) + ba(x,t)—(f)(x) + — a 2 (x,t)-^<f>(x) 

d 



+ 



(x + ya(x,t)) -4>{x) (y)ya(x,t)—(f>(x) 



R\{0} 



^(dy), 
(14) 
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where the indicator function of the set (—1,1). To get the second 

identity in (fT5j) . we have used the fact that 



E 



f r t+At d 

< J a(X s -,s)—ct>(X s )dB t 



X t = x } = 0, 



(15) 



E 



{ rt+At /• 
/ / [<P(X s .+ya(X s 
Jt J\y\<l 



N(ds,dy) 



X t = x \ = 0, 



(16) 



and 
E 



{rt+At r 
/ / [</>(X a -+ya(X a -))-<f>(X a -,8)] N(ds,dy) 
Jt J\y\>l 



X t = x 



t+At 



< (x + ycr(x, s)) — <fi(x)] v(dy)ds. 



(17) 



It J\y\>l 

Denoting u(x,t) = E ((ft(Xt) \ Xq = x), it can then be shown that [T] 

d 

—u(x,t) = Qu(x,t), 



(18) 



which is the well known backward Kolmogorove equation. Let p(y, t) be the 
probability density function for X t associated with the SDE ©. Then (fT8j) 
becomes 



d 

Oi(<t>(y)p(y,t)) dy= I Q4>(y)p(y,t)dy. 



Let Q* be the adjoint operator of Q. Therefore, 



(19) 



Q4>{y)p(y,t)dy = \ <j)(y)Q*p(y,t)dy 



(20) 



With relation (J20J) in mind, equation (|19p becomes 







<i>(y) ( Q- t p(v,t) - Q*p(y,t) ) dy = o. 



Since (f2T|) is true for any (p € Cq°(R), we get 

d 



(21) 



(22) 



which is the Fokker-Planck equation, the governing equation for the transi- 
tion probability density p for the stochastic dynamical system 1^ . 

Now we try to derive an expression for the adjoint operator Q* . It follows 
from (|14p that the operator Q can be written as 



Q = Ai + A 2 



(23) 



where A\ and A 2 are defined as 

A 1( p{x) = f( x ,t)^-(f>(x) + ba(x,t)^-<P{x) + ^a 2 {x,t)^(j){x), (24) 



and 

A 2 Mx) 



d 

4>(x + ya(x,t)) - <p(x) - I^ 1A )(y)ya(x,t)—<p(x 



R\{0} 



^(dy), 
(25) 



respectively. Note that A\ is expressed as 
A* lP (x,t) -- 



|- [p{x, t) (f(x, t) + ba{x, t))] + ( cj2 ( x ' 'M*. *))■ ( 26 ) 



We now find an expression for A\. Once A\ is obtained, the adjoint operator 
Q* can be expressed as 



Q* = A\ + A* 2 

and then the Fokker-Planck equation pt = Q*p will be obtained. 
By Taylor expansion, (f)(x + ya(x)) can be written as 

OO u 

4>{x + ya{x, t)) = <f>(x) + Y\ ak ( X ' ^dxk^' 



(27) 



fe=i 



(28) 



Substitute (f28|) into (fl~4"j) and ([25]) . respectively, we get 



+ 



fl\{0} 



,fc=i 



f(dy), 
(29) 



and 



Ao6(x) 



R\{0} 



7c=l 







u(dy). 



(30) 



It follows from (|30p that 

JO 

j42</>(x)p(x, i) dx 



i\{0} 



5^ — a fc (a;, t)-^$(x) - ^(-1,1) (2/)^ (x, *) 



K\{0} \K 



,fc=l 
00 „,fc 



(9 



(9^ 



Lk=l 

00 C „^ft 



K\{0} \K 



k=l 



dx 



v{&y)j p(x,t)dx 
p(x,t)dx S \ v{dy) 

x) dx I v{dy) 



,fc=i 



dx 



l\{0} 

(x, t)A^ p(x, t) dx, 



where A\, the adjoint of A2, is as follows 



v{dy) J 0(x) dx 
(31) 



A* 2 p(x,t) 



,fc=i 



(-y) k d k ( 

kl dx k 



a (x,t)p(x,t)j + I { _ 1A) (y)y—((r(x,t)p(x,t)) 



v(dy). 



(32) 



Note that in obtaining the third identity in f)31 [) . we have made use of 

p(x, t)a k (x, t)^<t>{x) dx = {-if J <t>{x)^ *)) dx > yk G N - 

R 

(33) 



Here N is the set of natural numbers. Given (|32 j) . {J26J) and (|27p . it follows 
from (I22p to get the desired Fokker-Planck equation. We summarize this 
result in the following theorem. 
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Theorem 1. Under the Assumption (Hi), the Fokker-Planck equation 
for the ltd SDE (0) is 

dv(x t) d d Id 2 

y = (f(x,t)p(x,t)) - b—(a(x,t)p(x,t)) + -A—{a 2 {x,t)p(x,t)) 



R\{0} 



E 

Lfc=l 



<a k (x, t)p(x, t)) + /(_!,!) (y)y— (a(x, t)p(x, t)) 



kl dx k 



(34) 



Remark 1. Two special cases are noted here. 

(i) When a(x,t) = 1, the Fokker-Planck equation |3^[ ) reduces to 

dj ^r = -it ^ t)p{x > t)} - t)} + w t} 



iJ\{0} 



p(x-y,t) - p(x,t) + I(-i,i)(y)y-^p(x,t) \ u(dy). 



(35) 



27ms is i/ie Fokker-Planck equation for an ltd SDE with an additive Levy 
process. 

(ii) When the Levy jump measure v satisfies f^ >1 \x\ k v(dx) < oo for 
k £ N, which is true if and only if~E(\Xt\ k ) < oo (see fl^), the Fokker-Planck 



equation \34\) can be rewritten as 



= -J^(/MMM)) + f>feJ^ (a k {x,t)p{x,t)) , (36) 



k=l 



where 



ci 



yv(dy), 



>i 



(37) 



and 



C2 



Cfc 



2 2 



y 2 v(dy), 



v(dy), k > 3. 



(38) 



(39) 



9 



3 Fokker-Planck equations for SDEs with Marcus 
integrals 

In this section, we derive the Fokker-Planck equation for Marcus SDE ([U]), 
in which Lj being a scalar Levy process with the triplet as (b,l,u). For 
simplicity, we make the following assumption. 

Assumption (H2) 
Assume that the noise intensity a(x, t) = <j(x) 7^ 0. 

Then ([6]) becomes 

dX t = f(X t ,t)dt + a(X t )odL t . (40) 

Define 

Y t = H(X t ) = / —du. (41) 
Jo °W 

The transform H in (|41|) is the Lamperti transform [2j. 
By the chain rule of the Marcus integral [1] , we have 



dY t = H\X t ) odX } 



H'(X t )f(X t ,t)dt+(H'(X t )a(X t )) odL t 
f(X t ,t 



f(H~ 1 (Y t ),t) 



°{H-\Y t )) 



dt + dL t 

dt + dL t . (42) 



Introducing 



f{Yut) ~ v(H-HY t ))' (43) 

equation (|4"S|) becomes 

dY t = f(Y t ,t)dt + dL t . (44) 
By Theorem [2] in the last section, the Forkker-Planck equation for Yt is 

|p(y, t) = ~ (Ry, t) P (y, tj) - b^ P (y, t) + —1*31, t) 

+ J (p(y-r)-p(y) + I(-i t i)(r)r-j^p(y,t)Jv(dr). (45) 

R\{0} 
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Let q(x,t) represents the probability density function of Xf, from the 
transform (|41|) and the Fokker-Planck equation (|45|) . we get the following 
result for the desired Fokker-Planck equation for SDE (|40p defined by Marcus 
integrals. 

Theorem 2. Under the Assumptions (Hi) and (H2), the Fokker-Planck 
equation for the Marcus SDE fl^0[ ) is 

d d d A d ( d 

^{x) — q(x,t) = -a(x) — (f{x,t)q(x,t)) - ba(x)— (a(x)q(x,t)) + TT "^)^ ( a ^g^ ( a ( x )<l( x ^)) 

+ J \^{H~ 1 (H{x) - r))q(H^ 1 (H(x) -r),t)- a(x)q(x,t)) +/(_ lj i)(r)m(ar)^ (a(x)q(x,t)f) u(dr). 

R\{0} 

(46) 



Remark 2. Consider a special case. When a(x) = 1, UW reduces to 



dq(x,t) d , , , , d , Ad 2 , . 

-— (f(x,t)q(x,t)) - b—q(x,t)) + -—5 q(x,t) 



dt dx ' ' cte ' 2 <9x 2 

+ 

R\{0} 

(47) 



U(x-y,t) -q{x,t) (y)y-j^q(x,t)\ u(dy), 



which is equivalent to \35\). the Fokker-Planck equation for a SDE with an 
additive Levy process. This is the consequence of the fact that for additive 
Levy noise, the Marcus integral and Ltd integral are the same. 



4 Examples 

In this section, we present the Fokker-Planck equations for some SDEs with 
special Levy processes, such as a Brownian motion together with a Poisson 
process, a compound Poisson process, and finally a a-stable Levy motion, 
respectively. 

For simplicity, in all these examples, we take f(x, t) = x, and cr(x, t) = x. 
Then the Ito SDE © and the Marcus SDE © become 

dX t = X t dt + X t dL t , (48) 
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and 

dX t = X t dt + X t odL t , (49) 

respectively. 

Example 1. When Lt is a standard Brownian motion together with 
a Poisson process with parameter A, the triplet (b, A, v) of Lt is (0, 1, A<5i), 
where 5\ is the Dirac measure defined as 

[l, ifleA 

Si(D) = \ (50) 
[o, ifl^D. 

By using (|34p . we get the Fokker-Planck equation for (|48|) 

fc=l 

(51) 

By using (06]), we get the Fokker-Planck equation for (|49|) 
<9<?(x,t) 



0i 



x— (xg(x, t)) H — x— ( x— (xp(x, t)) J + Ae Mxlq^e 1 |x|, i) — Xxq(x, t). 
ox 2 ox \ ox J 

(52) 



Example 2. When Lj is a compound Poisson process defined as 

N t 

L t = ^i, (53) 

i=l 

where Nt is a Poisson process with parameter A and (i = 1, 2, • • • ) are 
i.i.d. random numbers with standard normal distribution N(0, 1), the triplet 
(b, A, v) of L t is b = 0, A = 0, and z/(dx) = -^^e 2 dec. By using (|34"|) . we 
get the Fokker-Planck equation for (|48p 

^r 1 = -^(/(MMx,t))+Af;c^ , (54) 

fe=i 

where 

0, if A: is odd, 

c k = { (k€ N), (55) 

(k — 1)!!, if k is even. 



12 



71!! 



(56) 



and !! is the double factorial denned by 

n ■ (n — 2) ■ ■ ■ 1, if n is odd, 
n ■ (n — 2) • • • 2, if n is even. 

By using (I46p . we get the Fokker-Planck equation for (I49p 
x — -j — = — x— — (xg(x, i)) + A / (e~ r \x\q(e~ r \x\, t) — xq(x, t)) dr. 

ut OX 7R\{0} 

(57) 

Example 3. When Lj is a symmetric a-stable Levy motion with the 
triplet as b = 1, A = and u(dx) = ; where C Q is a constant (de- 

pending on a). Then by using (|34p . we get the Fokker-Planck equation for 



(— y) fc <9 fc 
A;! dx k 







a '(x,t)p(x,t)j + /(_i,i)(y)y^ (cr(x,t)p(a;,t)) 



(58) 



ll+O! ' 



By using (|46|) . we get the Fokker-Planck equation for (|49|) 



dq(x, t) 

+ xc a 



x iL ( xq ( x,t ^ 







\x\q(e r \x\,t) - xq(x,t)I { _ ltl) (r)r— (a(x,t)p(x,t)) 



dr 



\l+a " 



(59) 
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